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ABSTRACT 

The evolution of a quantum lattice gas automaton (LGA) for a single charged particle is 
invariant under multiplication of the wave function by a global phase. Requiring invari- 
ance under the corresponding local gauge transformations determines the rule for minimal 
coupling to an arbitrary external electromagnetic field. We develop the Aharonov-Bohm 
effect in the resulting model into a constant time algorithm to distinguish a one dimen- 
sional periodic lattice from one with boundaries; any classical deterministic LGA algorithm 
distinguishing these two spatial topologies would have expected running time on the order 
of the cardinality of the lattice. 
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Quantum lattice gas automata (QLGA) have been proposed as a possible architecture 
for solid state quantum computers since they require only an array of sites which can sup- 
port an extended (multi-)electron wave function [1]. The simplicity of such an architecture 
makes nanoscale fabrication plausible [2]. The main incentives for pursuing the program 
of quantum computation, however, are the quantum algorithms of Shor [3] and Grover 
[4], for example, which provide substantial improvements over deterministic or probabilis- 
tic algorithms. These quantum algorithms can be efficiently implemented with quantum 
gate arrays [5] , but at least in simple translations — analogous to deterministic billiard ball 
models for universal computation [6] — seem impractical to implement in QLGA models. 
Just as deterministic LGA efficiently simulate fluid flow in certain parameter regimes [7], 
QLGA seem best suited for simulation of quantum physical processes [1,8]. Single particle 
QLGA have been shown, in fact, to limit to the Schrodinger [9] and Dirac [1] equations 
under appropriate conditions. In this letter we show that in simulating a single quantum 
particle, QLGA efficiently perform a new and interesting computation — of spatial topology. 

We draw inspiration from the topological invariance of the spectral flow in a one 
parameter family of Dirac operators [10]. Spectral flow derives from the components of the 
group of gauge transformations, so after recalling the definition of one dimensional single 
particle QLGA with inhomogeneous potentials, we analyze which evolution rules are gauge 
equivalent. We find a difference between the gauge equivalence classes of rules on periodic 
lattices and lattices with boundary, and show how this may be exploited to distinguish 
between these two spatial topologies. We conclude with a discussion of the complexity of 
this quantum computation and point out directions for further investigation. 

A one particle QLGA is a discrete time model for a quantum particle moving in an 
array (lattice) of sites. Here we will consider only finite one dimensional lattices L, with 
or without periodic boundary conditions. In these cases we need only allow two velocities 
{±1} in order to construct a model which limits to the Schrodinger [9] or Dirac [1] equation 
as the timestep and lattice spacing scale to zero appropriately. The amplitudes for the 
particle to be (left, right) moving at a lattice site x G L combine into a two component 
complex vector tp(t, x) := (-i/>-i(t, x), ip+i(t, x)) which evolves as 



Here the weights wp(x) £ M 2 (C) are 2x2 complex matrices constrained by the requirement 
that the global evolution matrix 



ip(t + 1, x) = W-i(x)i()(t, x — 1) + w+i(x)ip(t, x + 1). 



\ 









U := 



w-i(x) 





W-i(x + 1) 









w +1 (x + 1) 



V 



be unitary. 
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For periodic boundary conditions the top and bottom rows (of blocks) of U wrap 
around and completely homogeneous solutions are possible. We showed in [1] that the 
most general parity invariant homogeneous solutions, up to unitary equivalence and an 
overall unobservable phase, form a one parameter family 

, N / i sin 9 \ , N / cos 9 \ 

w-i(x) = w-i := „ = tu+i := ../) n , 

v y ^0 cost' y \ v «sin6' Oy 

where 6* scales to the mass in the Dirac equation limit. 

Without periodic boundary conditions, the local rule must change at the boundary. 
We show in [11] that there is a one parameter family of (Type II) boundary conditions 
which give the global evolution matrix the form: 



U 



W-i 



w+i 


W-i 



W+i 





(2) 



at the left boundary. Here 



w+i ■= 




for any ( G 1R and the right boundary condition is characterized by «J_i obtained by 
parity transformation, although with an independent phase. For simplicity we set both 
phase angles to zero in the following discussion; this will not affect the relevant features of 
our results. 

The global evolution matrix U acts by left multiplication \I>(£ + 1) = U^(t) on the 
wave function *(£) : L x {±1} — ► C. Multiplication of the wave function by an overall 
phase e~ lOL leaves the quantum state invariant since the phase cancels in any observable 
(\E r |(9|\E r ). Requiring invariance of the evolution under the local gauge transformations 
e -ia(t,x) corresponding to this global symmetry necessitates modification of U: Setting 
yj'(t, x) := e~ ia ( t,x ^yj(t, x), which we write as 

V'(t) :=D[e-* a(t ' x) ]*(t), (3a) 
we achieve gauge invariance with 

U'{t) := D[e- ia{t+1 ^]UD[e ia ^ x) ] (36) 

since 

U'(t)y'(t) = D[e- ia(t+1 ' x) }UD[e ia{t ' x) ]D[e- ia(t ' x) ]^(t) 
= D[e- ia{t+1 ' x) ]^(t + l) 
= *'(£+!). 
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Expanding (36) we find that the superdiagonal blocks of U'(t) are 

U'(t) x - ljX = e -i<x(t+i,x-i)+ia(t,x) w+u ( 4a ) 
while the subdiagonal blocks are 

U'(t) x , x -! = e -Mt+l,x)+ia(t,x-l) w _^ ( 4fe ) 

For a periodic lattice, x runs over all the labels {0, . . . , \L\ — 1} of the lattice sites and the 
indices are interpreted mod \L\. For a lattice with boundaries, x E {1, . . . ,\L\ — 1} since 
there are no blocks in the evolution matrix containing nonzero amplitudes for transitions 
between x = and x = \L\ — 1. Multiplying and dividing each block in column x (with 
these same conventions) of U'(t) by e ia ( t + 1 ' x ) j the expressions in (4) can be rewritten as 

U'{t) x ^ x = e -i*ta(t+l,x)+iA x a(t+l,x) w+i ^ 

and 

U\t) x , x - X = e -^Mt+hx-i)-iA x a(t+i,x) w _^ 

where A is the difference operator so A t a(t + 1, x) := a(t + 1, x) — a(t, x) and similarly 
for A x . In column x both the super- and subdiagonal blocks now contain the phase factor 
e -iA t a(£+i,x) R eca vi that when we included an inhomogeneous potential in a QLGA [12] 
it also had the effect of multiplying the blocks in column x by a phase factor, e~ % ^ x > * 
Thus the gauge transformation (3) transforms a potential (f)(t, x) to 

(f)'(t,x) = </>(t,x) + A t a(t+ l,x). (6) 

The A x a phases in the gauge transformed evolution matrix U'(t) remind us that the 
electromagnetic potential has a second component — the vector potential A. The conjugate 
A x a phases of the antidiagonal blocks in (5) indicate that the vector potential enters the 
evolution matrix in the same way. That is, in the presence of an electromagnetic potential 
(4>(t, x), A(t, x)), the super- and subdiagonal blocks in the evolution matrix are 

U^ A \t) x - ltX = e -i*(t,*)+iMt,x) w+1 ( ?a ) 

and 

U^ A \t) x ^ = e -i*(t,*)-iMt,x) w _ u (76) 

respectively. Notice that for a lattice with boundaries, A(t, 0) = A(t, \L\) does not en- 
ter into these expressions for the global evolution, although it does for a periodic lattice. 
jj((t>,A) degcribeg the discrete version of the usual minimal coupling of a Dirac particle to an 



Column and row are interchangeable here, corresponding to multiplying by the local potential phase 
factors before or after the advection part of each timestep, respectively. The algorithm described 
in [12] multiplies before, while those in [9] and [13] do so afterwards. To obtain the expressions 
corresponding to those in (5), but with a common A t a phase in each row, we can multiply and 

divide each block in row x of U'(t) by e ta(t ' x) . 
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external (classical) electromagnetic field [14]. The conjugate phases from the vector poten- 
tial are the consequence of the velocity dependence of this coupling; the vector potential 
itself transforms to 

A'(t,x) = A(t,x) + A x a(t + l,x) (8) 

under the gauge transformation (3), analogously to the transformation (6) of the scalar 
potential. 

The consequences of this gauge freedom in the vector potential differ for periodic 
lattices and those with boundaries: Given any vector potential A(t, x), consider the gauge 
transformation defined by 

X 

a(t + i,x) = -J2Mt,y)- (9) 

With this definition, A x a(t + 1, x) = —A(t, x) for 1 < x < \L\ — 1 so that using (9) in (8) 
gives a gauge transformed A'(t, x) which vanishes except at x = 0. Since the global evo- 
lution on a lattice with boundaries does not depend on A'(t, 0), (9) gauge transforms the 
vector potential away completely — there is exactly one gauge equivalence class of vector 
potentials. On a periodic lattice, however, the global evolution depends on A'(t, 0) (mod 
27r), so there is a one (periodic) parameter family of gauge equivalence classes. Alterna- 
tively, the gauge transformation (9) demonstrates that the only gauge invariant quantity 
which can be constructed from the vector potential is (a function of) 

\L\-1 

5 := expz A (t, x ), ( 10 ) 

upon which the evolution depends only for 
a periodic lattice, in which case (10) is the 
discrete analogue of a Wilson loop vari- 
able/ho lo no my, i.e., a winding number. 

Such a gauge invariant quantity man- 
ifests itself in the spectrum of the evolu- 
tion operator: As the phase 5 varies, so 
does the spectrum of — if the lattice 

is periodic, but not if it has boundaries. 
Figure 1 shows an example of this spectral 
flow [10] for a periodic lattice with N = 16 
and 6 = n/6. Here 4>(t : x) has been set to 
by a gauge transformation, and 5 varies 
from to 2-7T. As 5 increases, the fre- 
quencies (energies) of positive frequency 
right /left moving plane waves, i.e., the 
eigenvectors of U^' A \ increase/decrease, 
respectively. When the 'mass' 9 = 0, half 



CO 




Figure 1. Spectral flow in a periodic QLGA with 
\l\ = 16 and = ir/6. The scalar potential 
has been taken to vanish and 5, the gauge invari- 
ant degree of freedom in the vector potential runs 
from to 2ir. 
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the eigenvalues shift up by 1 and half shift down by 1 as 6 varies from to 2n; the spectral 
flow (the number of eigenvalues crossing each frequency level) is 2. 

To observe such a spectral flow, imagine preparing a quantum particle on a periodic 
lattice in a positive frequency, right moving eigenstate of U^ 0,0 \ but then applying an 
external vector potential so that the evolution is by U^ '^ instead, where A(t, x) = A. 
A massless particle will remain in the original eigenstate, but a frequency (energy) mea- 
surement at any subsequent timestep will return the original value plus A = S/\L\. For a 
massive (i.e., 9 ^ 0) particle, the initial state is a superposition of the positive and negative 
frequency eigenstates of U^°' A ^ with the original wave number, but for small masses the 
positive frequency still dominates so again a frequency (energy) measurement will return 
the original value plus an amount on the order of A. 

Our goal, however, is to distinguish between a periodic lattice and one with bound- 
aries. Since the single particle eigenstates differ for these two situations even if \L\ is the 
same [13,11], without knowing the spatial topology a priori we cannot prepare a quantum 
particle in an eigenstate. But we can prepare a positive frequency right moving wave packet 
with bounded support — identically on either lattice. Subsequent frequency measurements 
will return a distribution of frequencies concentrated around the expectation value. This 
distribution is the same for either lattice, if the evolution is by U^°'°\ But if the evolution 
is by U ( - 0,A ' ) with A > 0, say, the observed frequency distribution will be shifted to larger 
values on the periodic lattice. 

Since some finite, fixed, number of measurements suffices to distinguish these distribu- 
tions with probability 1 — e for a given initial wave packet and external vector potential, the 
computation takes only constant time, independent of the size of the lattice \L\. In fact, 
for larger lattices the initial wave packet can be broader, i.e., more concentrated around 
its expected energy, and hence fewer measurements are needed to distinguish the original 
distribution from the one shifted by the vector potential in case the lattice is periodic. In 
contrast, suppose we try to distinguish a periodic lattice from one with boundaries using 
a deterministic LGA. Since no superpositions of states are possible, all that we can do is 
to start a single particle off in one direction and see if it ever changes velocity by reflecting 
from a boundary. To do so the particle must travel for 0(|L|) timesteps, on the average, 
so the QLGA algorithm provides a even greater improvement over the classical algorithm 
than does Grover's algorithm for searching [4]. 

In conclusion, we remark that our algorithm exploits the Aharonov-Bohm effect [15] 
which is more usually discussed in two dimensions — we are imagining creating the vector 
potential by applying a magnetic field which threads a (possibly incomplete, in the non- 
periodic case) ring of lattice sites. It would be interesting, and potentially useful for pattern 
recognition [16], to formulate this algorithm for a two dimensional QLGA. Preparation of 
a localized wave packet, and measurement of its energy, each in constant time, is then 
plausible on physical grounds, provided that the lattice lies within some fixed area. For 
the more realistic situation of a fixed density of lattice sites (and hence an increasing spatial 
area with increasing lattice size) additional analysis is required. As the example of unitary 
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transformations and measurements on the Rydberg states of an atom illustrates, careful 
attention to the details of the physical implementation of a quantum algorithm is required 
to correctly quantify its computational complexity [17]. Although the evolution of the 
QLGA can be simulated in the standard poly-local model for quantum computing [18], 
issues of the time required for state preparation and measurement, and the adiabaticity 
required when turning on the magnetic field, all must be considered before we can claim 
to understand how this computation scales physically. 
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